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, Abstract. We introduce a model for self-organized dynamics which, we argue, addresses 

several drawbacks of the celebrated Cucker-Smale (C-S) model. The proposed model does 
not only take into account the distance between agents, but instead, the influence between 
agents is scaled in term of their relative distance. Consequently, our model does not 
involve any explicit dependence on the number of agents; only their geometry in phase 
space is taken into account. The use of relative distances destroys the symmetry property 
of the original C-S model, which was the key for the various recent studies of C-S flocking 
behavior. To this end, we introduce here a new framework to analyze the phenomenon 
, of flocking for a rather general class of dynamical systems, which covers systems with 

non-symmetnc influence matrices. In particular, we analyze the flocking behavior of the 
proposed model as well as other strongly asymmetric models with "leaders". 
, The methodology presented in this paper, based on the notion of active sets, carries 

' over from the particle to kinetic and hydrodynamic descriptions. In particular, we discuss 

the hydrodynamic formulation of our proposed model, and prove its unconditional flocking 
for slowly decaying influence functions. 
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1. Introduction 

The modeling of self-organized systems such as a flock of birds, a swarm of bacteria or a 
school of fish, [H m O [121 IHl EOl [ZD [26], has brought new mathematical challenges. One 
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of the many questions addressed concerns the emergent behavior in these systems and in 
particular, the emergence of "flocking behavior". Many models have been introduced to 
appraise the emergent behavior of self-organized systems [11 [3l [71 [131 [13 [221 [23 [2Z] • The 
starting point for our discussion is the pioneering work of Cucker-Smale, [HI [9], which led to 
many subsequent studies [3l [6l [HI [151 [161 [22] . The C-S model describes how agents interact 
in order to align with their neighbors. It relies on a simple rule which goes back to [22]: the 
closer two individuals are, the more they tend to align with each other (long range cohesion 
and short range repulsion are ignored). The motion of each agent "i" is described by two 
quantities: its position, Xj G W^, and its velocity, Vj E W^. The evolution of each agent is 
then governed by the following dynamical system, 

(1.1a) — =v„ -^ = ^2^M-^--^)- 

i=i 

Here, a is a positive constant and (pij quantifies the pairwise influence of agent "j" on the 
alignment of agent "i" , as a function of their distance, 

(1.1b) := 0(|xj - 

The so-called influence function, is a strictly positive decreasing function which, by 
rescaling a if necessary, is normalized so that i;^(0) = 1. A prototype example for such an 
influence function is given by (p{r) = (1 + r)~^, s > 0. Observe that the C-S model (jl.ip is 
symmetric in the sense that the coefficients matrix (pij is, namely, agents "i" and "j" have 
the same influence on the alignment of each other, 

(1.2) (f)ij = 

The symmetry in the C-S model is the cornerstone for studying the long time behavior of 
(jl.ip . Indeed, symmetry implies that the total momentum in the C-S model is conserved, 

(l-3a) ^(^E^*(*))=0 ^ v(t):=l5^v,(0 = v(0). 

\ i=l / i=l 

Moreover, the symmetry of (jl.2p implies that the C-S system is dissipative, 

(1.3b) |v,- — vP = ^ (t>7iWi — VoP < — mmSjiit) x — |v,- — vP. 

t i,j t 

Consequently, (|1.3p yields the large time behavior, Xi(t) k, vt, and hence miuj^ (/)jj(t) > 
(/)(|v|t). This, in turn, implies that the C-S dynamics converges to the bulk mean velocity, 

(1.4) v,(i)*^v(0), 

provided the long-range influence between agents, </>(|xj — Xj|), decays sufficiently slow in 
the sense that (/>(•) has a diverging tail, 

/oo 
0(r) dr = oo. 

We conclude that the C-S model with a slowly decaying influence function (jl.Sp . has an un- 
conditional convergence to a so-called flocking dynamics, in the sense that (i) the diameter, 
maxjj- |xj(t) — Xj{t)\, remains uniformly bounded, thus defining the domain of the "flock"; 
and (ii) all agents of this flock will approach the same velocity — the emerging "flocking 
velocity" . 
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Definition 1.1. |151 p. 416] Let {xj(t),Vj(t)}j=i^ be a given particle system, and let 
dx{t) and dv{t) denote its diameters in position and velocity phase spaces, 

(1.6a) dx{t) = max |xj(t) — Xj(t)|, 

(1.6b) dY{t) = max |vj (t) — Vj(t)|. 

The system {xj(t), Vj(t)}j=i^...^Ar is said to converge to a flock if the following two conditions 
hold, uniformly in N, 

(1.7) sup dx{t) < +00 and lim dv{t) = 0. 

t>0 t-^+oo 

Remark 1.2. One can distinguish between two types of flocking behaviors. When (|1.7|) 
holds for all initial data, {xi(0), Vj(0)}j=i^...^Ar, it is referred to as unconditional flocking, 
e.g., [51 [HI [T ^ I15[ \2'6\. In contrast, conditional flocking occurs when (|1.7|) is limited to a 
certain class of initial configurations. 

The flocking behavior of the C-S model derived in [15] was based on the £^-based argu- 
ments outlined in (jl.3p . Other approaches, based on spectral analysis, i^- and ^°°-based 
estimates were used in [6l [H [T3j to derive C-S flocking with a (refined version of) slowly 
decaying influence function (jl.5p . Though the derivations are different, they all require the 
symmetry of the C-S influence matrix, 

Despite the elegance of the results regarding its flocking behavior, the description of self- 
organized dynamics by the C-S model suffers from several drawbacks. We mention in this 
context the normalization of C-S model in (ILlah by the total number of agents, A^, which is 
shown, in section 12.11 below, to be inadequate for far-from-equilibrium scenarios. The first 
main objective of this work is to introduce a new model for self-organized dynamics which, 
we argue, will address several drawbacks of the C-S model. Indeed, the model introduced in 
section [22] below, does not just take into account the distance between agents, but instead, 
the influence two agents exert on each other is scaled in term of their relative distances. As 
a consequence, the proposed model does not involve any explicit dependence on the number 
of agents — just their geometry in phase space is taken into account. It lacks, however, 
the symmetry property of the original C-S model, ()1.2p . This brings us to the second main 
objective of this work: in section [3] we develop a new framework to analyze the phenomenon 
of flocking for a rather general class of dynamical systems of the form, 

dxj dvj ■r-^ 

^ = v„ ^ = « 2^ - v.), aij > 0> 2^ < 1, 

which allows for non-symmetric influence matrices, a^j ^ aji. Here we utilize the concept 
of active sets, which enables us to define the notion of a neighborhood of an agent; this 
quantifies the "neighboring" agents in terms of their level of influence, rather than the 
usual Euclidean distance. The cornerstone of our study of flocking behavior, presented in 
section 13.11 is based on a key algebraic lemma, interesting for its own sake, which bounds 
the maximal action of antisymmetric matrices on active sets. Accordingly, the main result 
summarized in theorem 13. 4[ quantifies the dynamics of the diameters, dx{t) and dv{t), in 
terms of the global active set associated with the model. We conclude, in section [J that 
the dynamics of our proposed model will experience unconditional fiocking provided the 
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influence function (p decays sufficiently slowly such that, 

/oo 
(fP'{r)dr = oo. 

This is slightly more restrictive than the condition for ffocking in the symmetric case of C-S 
model, (jl.Sp . Another fundamental difference between the ffocking behavior of these two 
models is pointed out in remark 1412) below: unlike the C-S ffocking to the initial bulk velocity 
v(0) in ()1.4p . the asymptotic ffocking velocity of our proposed model is not necessarily 
encoded in the initial configuration as an invariant of the dynamics, but it is emerging 
through the ffocking dynamics of the model. 

The methodology developed in this work is not limited to the new model, whose flocking 
behavior is analyzed in section [^TTl In section W?2\ we use the concept of active sets to study 
the ffocking behavior of models with a "leader". Such models are strongly asymmetric, 
since they assume that some individuals are more influential than the others. 

Finally, in section [5] and, respectively, section [6l we pass from the particle to kinetic and, 
respectively, hydrodynamic descriptions of the proposed model. The latter amounts to the 
usual statements of conservation of mass, p, and balance of momentum, pu, 

(1.9a) + Vx • (pu) = 

(1.9b) at(pu) + Vx(pu0u) = ap(^^-u^ , (w)(x) :=^(/.(|x-y|)w(y)p(y)(iy. 

We extend our methodology of active sets to study the ffocking behavior in these contexts 
of mesoscopic and macroscopic scales. In particular, we prove the unconditional flocking 
behavior of (II. 9p with a slowly decaying influence function, such that (jl.Sp holds, 

sup |u(t,x) - u(t,y)| *^ 0. 

x,yGSiipp(p(-,t)) 



2. A MODEL FOR SELF-ORGANIZED DYNAMICS 

In this section, we introduce the new model that will be the core of this work. This model 
is motivated by some drawbacks of the C-S model. 

2.1. Drawbacks of the C-S model. Originally, the C-S model was introduced in [8] to 
model a ffnite number of agents. The normalization pre- factor 1/A^ in (jl.lap was added 
later in Ha and Tadmor, [15], in order to study the "mean-fleld" limit as the number of 
agents becomes very large. This modiffcation, however, has a drawback in the modeling: 
the motion of an agent is modifled by the total number of agents even if its dynamics is only 
influenced by essentially a few nearby agents. To better explain this problem, we sketch a 
particular scenario shown in flgure[TJ Assume that there is a small group of A^i agents, Gi, 
at a large distance from a large group of N2 agents, G2; by assumption, we have A^i « N2. 
If the distance between the two groups is large enough, we have, 

(2.1a) if i € Gi and j G G2. 

In this situation, the C-S dynamics of every agent "z" in group Gi reads, 
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Therefore, since there are only A''i "essentiahy" active neighbors of "i" , yet we average over 
the much larger set of A''! + ^ -^i agents, we would have dvi/dt ?s 0. Thus, the presence 
of a large group of agents G2 in the horizon of Gi, will almost halt the dynamics of Gi. 




Figure 1. A small group of birds Gi at a large distance from a larger group 
G2 (|2.1ap . Due to the normalization in the C-S model (jl.lap . the group 
Gi will almost stop interacting. 



2.2. A model with non-homogeneous phase space. We propose the following dynam- 
ical system to describe the motion of agents {xj(t), Vj(t)}^-^, 

(2-2) ^ = ^ = ^n" , "^^J ~ ^i)' "^^J = '^(I^'J " ^^\^- 

2-^k=\ rik j=l 

Here, a is a positive constant and is the influence function. The main feature here is 
that the influence agent "j" has on the alignment of agent "i", is weighted by the total 
influence, 'Yl,k=i't'iki exerted on agent "i". 

In the case where all agents are clustered around the same distance, i.e., w 0o, then 
the model (|2.2p amounts to C-S dynamics. 



avj a 



dt N(l)o . 

But unlike the C-S model, the space modeled by (j2.2p need not be homogeneous. In 
particular, it better captures strongly non-homogeneous scenarios such as those depicted in 
I2.lt the motion of an agent "i" in the smaller group Gi will be, to a good approximation, 
dominated by the agents in group Gi, 

dvi « i, / X 

Here, (po is the coefficient of interaction inside the nearby group Gi, i.e., (pij ~ (po for 
i, j € Gi, whereas the agents in the "remote" group G2, will only have a negligible influence. 
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The normalization of pairwise interaction between agents in terms of relative influence 
has the consequence of loss of symmetry: the model (I2.2p can be written as, 

^=v., — = a>^a.,(v,-v,), 

i=i 

where the coefficients a^j, given by, 

22k=in\^k-^i\) 

lack the symmetry property, Uij ^ aji. Two more examples of models with asymmetric 
influence matrices will be discussed below. The flocking behavior of a model with leaders, in 
which agents follow one or more "influential" agents and hence lack symmetry, is analyzed 
in section H] below. In section [7] we introduce a model with vision in which agents are 
aligned with those agents ahead of them, as another prototypical example for self-organized 
dynamics which lacks symmetry, and we comment on the difficulties in its flocking analysis. 
Tools for studying flocking behavior of such asymmetric models are outlined in the next 
section. 

3. New tools to study flocking 

We want to study the long time behavior of the proposed model ()2.2D . The lack of 
symmetry, however, breaks down the nice properties of conservation of momentum, (jl.3ap . 
and energy dissipation, (|1.3b|) . we had with the C-S model. The main tool for studying the 
C-S flocking was the variance, |vj — in either one of its ^^-versions, p = 1,2 or 

p = oo. But since the momentum is not conserved in the proposed model (j2.2p . the variance 
is no longer a useful quantity to look at; indeed, it is not even a priori clear what should be 
the "bulk" velocity, v, to measure such a variance. 

In this section, we discuss the tools to study the flocking behavior for a rather general 
class of dynamical systems of the form, 

(3.1a) ^ = Vj, = - Vj), Uij > 0. 

Here, a is a positive constant, and aij > quantifies the pairwise influence of agent "j" on 
the alignment of agent "i", through possible dependence on the state variables, {x^, v/jj^. 
By rescaling a if necessary, we may assume without loss of generality that the 
normalized so that 

(3.1b) X^^^i^l- 
Setting an := 1 — ^iji '^^ rewrite ()3.ip in the form 
(3.2a) — =v„ — = a(v,-v,), 

where the average velocity, Vj, is given by a convex combination of the velocities surrounding 
agent "z", 

N 

(3.2b) Vi:=^aijVj, y^Qij = 1- 
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We should emphasize that there is no requirement of symmetry, allowing aij ^ aji. This 
setup includes, in particular, the model for self-organized dynamics proposed in (j2.2p . with 
asymmetric coefficients a^j = (kij/Ylk ^ik- 

In order to study the flocking behavior of (j3.ip , we quantify in section 13.21 the decay of 
the diameter, dy[t) using the notion of active sets. The relevance of this concept of active 
sets is motivated by a key lemma on the maximal action of antisymmetric matrices outlined 
in section 13.11 This, in turn, leads to the main estimate of theorem 13. 4^ which governs the 
evolution of dxif) and dy{t). 

3.1. Maximal action of antisymmetric matrices. We begin our discussion with the 
following key lemma. 

Lemma 3.1. Let S he an antisymmetric matrix, Sij = —Sji with maximal entry \Sij\ < M . 
Let u,w ^ be two given real vectors with positive entries, Ui,Wi > 0, and let U, W denote 
their respective sum, U = Ui and W = Wj. Fix 9 > and let \{9) denote the number 
of "active entries" of u and w at level 9, in the sense that, 

X{9) = \A{9)\, A{9) := {j I Uj > 9U and wj > 9W } . 

Then, for every 9 > 0, we have 

(3.3) \{Su,w)\ < MUW {I- X'^{9)9^) . 

Remark 3.2. Lemma \3.1\ tells us that the maximal action of S on u,w, does not exceed 

\{Su,w)\ < MU W mm {l-X^ (9) 9^). 



which improves the obvious upper-bound, \{Su,w)\ < MUW. 
Proof. Using the antisymmetry of S, we find 

{Su, VU) = ^ SijUiWj = ^ X] {'^i'^j - UjWi) , 

and since S is bounded by M, we obtain the inequality, 

\{Su,w)\ < ~ ""j^il- 

The identity, |a — 6| = a + b — 2 min(a, b) for a, 6 > 0, implies 

\{Su, w)\ < — (^UiWj + UjWi — 2mm{uiWj, UjWi} 

(3.4) = m(uW -'^mm{uiUJj,UjWi}y 

id 

By assumption, there are at least A(^) active entries at level 9 which satisfy both inequalities, 

keA{9): Uk> 9U and Wk> 9 W. 
Therefore, by restricting the sum in (j3.4p only to the pairs of these active entries we find 

\{Su,w)\ < m(uW - mm{uiWj,UjWi}^ < M {UW - X^{9) ■9U ■9W) , 

and the desired inequality ()3.3p follows. □ 
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3.2. Active sets and the decay of diameters. The concept of an active set aims to 
determine a neighborhood of one or more agents in (j3.ip based on the so-called influence 
matrix, {uij}, rather than the usual Euclidean distance. The following definition, which 
applies to arbitrary matrices, is formulated in the language of influence matrices. 

Definition 3.3 (Active sets). Let {aij} be a normalized influence matrix, Oij > 0, Ylj ^ij — 
1. Fix > 0. The active set, A.p{0), is the set of agents which influence "p" more than 9, 

(3.5) Ap{e) := {j I flp, > 6}. 

The global active set, A{6), is the intersection of all the active sets at that level, 

(3.6) Ai9) = f]Ap{e). 

p 

This notion of active set, Ap{6), defines a "neighborhood" for agent "p", and can be 
generalized to more than just one agent. For example, 

(3.7) Ap,{e):=Ap{e)nAg{e), 

is the set of all agents whose influence on both, "p" and "(;", is larger than 6, see figure [2j 

The number of agents in an active set Ax{6) is denoted by Xx{9), e.g. \pq{0) = \Apq{6)\. 
The numbers {\q{0)}pq are difficult to compute for general 0's: one needs to count the 
number of pairs of agents in the underlying graph Q, which stay connected above level 6, 



(3.8) 



1 if agent "i" is influenced by "j" : aij > 0, 
otherwise. 



One simple case we can count, however, occurs when 9 takes the minimal value, 9 = 
minjj Uij. Then, the active sets Ap{9) includes all the agents, Ap{9)e=minij aij = {1) • • • j 
and since this applies for every "p", then Apq{9) and the global active set, A{9), include all 
agents, 

(3-9) m\0=u^na,, = ^- 



AiW Am(^) 




Figure 2. An illustration of active sets. Here, Ai{9) = {1, 2, 3} and 
A4(0) = {2, 3, 4}. The pairwise active set, A14 = Ai{9) n A4(0), consists 
of agents "2" and "3" . 

Armed with the notion of active set and with the key lemma 13.11 on maximal action 
of antisymmetric matrices, we can now state our main result, measuring the decay of the 
diameters dx{t) and dyit) in the dynamical system (|3.2p . 
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Theorem 3.4. Let {xj(i), Vj(t)}j be a solution of the dynamical system (j3.2p . Fix an 
arbitrary 6 > and let Xpq{6) be the number of agents in the active sets, Apg{6), associated 
with the influence matrix of ()3.2p . Then the diameters of this solution, dx(t) and dv{t), 
satisfy, 

(3.10a) l'^^^*) - '^^^^^ 

(3.10b) ^dv{t) < -a min Xl JO) e"^ dv{t). 



dt 



pq 



Since A{6) C Ap^q{9) then Xpq{9) > X{d) and (I3.10bp yields the following global version 
of the theorem above. 

Theorem 3.5. Fix an arbitrary 9 > and let X{6) be the number of agents in the global 
active set, A{9), associated with (j3.2p . Then the diameters of its solution, dx{t) and dy{t), 
satisfy, 

(3.11a) Jt"^^^^^ - '^^^^^ 

(3.11b) ^dvit) < -aX'^(9)9^dv(t). 

dt 



Proof of theorem \3.4\ We fix our attention to two trajectories Xp(t) and Xg(i), where p and 
q will be determined later. Their relative distance satisfies, 

d 



Xp - Xgp = 2(xp-Xg, Vp-Vg) < 2jXp-Xg||Vp-V, 



dt 



which implies, 

^|Xp(t)-X,(t)| <dy(t). 

Thus, (|3.10ap holds. Next, we turn to study the corresponding relative distance in velocity 
phase space, 

(3.12) ^'^^"'^'''^ ^ 2a (vp-Vg , Vp-Vg) 

= 2a(vp-Vg , Vp-Vq) - 2a|vp-v<yp; 

recall that Vp and Vp are the average velocities defined in (I3.1bp . Given that Yle'^ke = Ij 
the difference of these averages is given by. 



« J « 

Inserting this into (j3.12p . we find, 

(3.13) ^ |vp - Vgp = 2a I ^apjagi(vp-Vg, Vj-Vj) - |vp - Vg| 

To upper-bound the first quantity on the right, we use the lemma [3TT] with Ui = Opi, wi - 
and the antisymmetric matrix Si^j = (vp — Vg, Vj— Vj): since \Sij\ < dy, U = X^jWi 
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and W = Wi = 1, we deduce, 

I ^apiOgi(Vp-Vg, Vj-Vi) < dy{l - Apq(6')6'^). 

Here, \pq{9) is the number of agents in the active set Apg(0), 

-^P<?(^) = l{j I "pj > ^ and a^j > 
Therefore, the relative velocity Vp — Vg in (|3.13p satisfies, 

^|vp-v,|2<2a(4(l-A2^(^)0 - |vp-v,| 

In particular, if we choose p and q such that |vp(t) — Vg(t)| = ci\/(t), the last inequality 
reads, 

(3.14) 4- dlit) < -2a min XlJe) 9^ dl(t). 

at pq 

and the inequality ()3.10bp follows. □ 

Remark 3.6. Equation ()3.10bp tells us that the diameter in velocity phase space, dv{t), is 
decreasing in time. In fact, an even stronger statement holds, namely, if we let Vl{t) denote 
the convex hull of the velocities, := Conw ({vi(t)}j=i..jv), then is decreasing in 
time in the sense of set inclusion, 

(3.15) ^{ti) D n{t2) ifti < t2. 

Indeed, by convexity, Vj G Q{t) for any i, and consequently, if Vi is at the frontier ofQ, 
then the vector (vj — Vj) points to the interior of Q at Vj, see figure\^ More precisely, if n 
denotes the outward-pointing normal to 0, at Vi, then Vj • n = (vj — Vj) • n < Therefore, 
the frontier of Q{t) is a 'fence" [18] for the vectors Vj(t) and (j3.15p follows. 




Figure 3. At the frontier of the convex hull il, the vector (vj — Vj) points 
to the interior of Thus, for any outward-pointing normal vector n at Vj, 
we have: Vj • n = (vj — Vj) • n < 0. 

The bound of dv{t) implies that the spatial diameter of the flock, dx{t) grows at most 
linearly in time. Indeed, for agents "p" and "q" which realize the maximal distance, dx{t) = 
|xp(t) — Xq(t)|, we have 

j/xit) < |Vp(i) -v,(i)| < dvit), 
and hence dx{t) < dx{^) + dy(0)t. 
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Theorem 13.41 and 13.51 will be used to prove the flocking behavior of general systems of the 
type (j3.2p . The key point will be to make the judicious choice for the level 6 = 6{dx{t)), to 
enforce the convergence dy(t) — >■ through the inequalities p.lOp . p. lip . In this context 
we are led to consider dynamical inequalities of the form, 

(3.16a) ^"^^(^^ - 

(3.16b) ^'^^(^^ - -(^i^idx{t))dv{t). 

The long time behavior of such systems is dictated by the properties of i/'(') > 0. 

Lemma 3.7. Consider the diameters dx{t),dy{t) governed by the inequalities (|3.16p . where 
ip{-) is a positive function such that, 



(3.17a) dy(0) < / '>p{r)dr. 

Then the underlying dynamical system convergences to a flock in the sense that (jl.7p holds, 

supdxit) < +00 and lim dv{t) = 0. 

In particular, if ip{-) has a diverging tail, 

/oo 
ip{r)dr = oo, 

then there is unconditional flocking. 

Proof. We apply the energy method introduced by Ha and Liu [H]. Consider the "energy 
functional", £ = £{t), 

rdx{t) 

(3.18) £{dx,dv){t) :=dv{t) + a i;{s)ds. 



The energy £ is decreasing along the trajectory {dx,dv), 

^£{dx,dv) = dv + aip{dx)dx < —a'ip{dx)dv + aTp{dx)dv = 0, 
and we deduce that, 

(3.19) dv{t) - dv{0) < -a ip{s)ds. 

JdxiO) 

By our assumption (|3.17ap . there exists > (independent of t), such that, 

(3.20) dv{0) = a [ ij{s)ds, 



dxiO) 



and the inequality (|3.19p now reads. 



fd, rdx{t) rd, 

dv{t)<a I ■ilj{s)ds — a / 'ip{s)ds = a / il){s)ds. 



ldx{0) JdxiO) Jdx(t) 

Since dv{t) > 0, we conclude that we have a flock with a uniformly hounded diameter, 
(3.21) dx{t)<d^ forallt>0, 
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thus improving the hnear growth noted in remark [3^ The uniform bound on dx{t) in (j3.2ip 
impHes that the velocity phase space of this flock shrinks as the diameter dv{t) converges 
to zero. Indeed, the inequality ()4.5bp yields, 

-i-dvit) < —aijj^, ■ dv(t), ■0* := ™in i^i^) > 0, 
dt 0<r<d, 

and Gronwall's inequality proves that dv{t) converges exponentially fast to zero. □ 

4. Flocking for the proposed model 

In this section we prove that the model (j2.2p converges to a flock under the assumption 
that the pairwise influence, 4>{\xj — Xj|), decays slowly enough so that has a non 

square-integrable tail, (jl.Sp . f°° (fp'{r)dr = oo. In section HT2l we show that the same result 
carries over the dynamics of strongly asymmetric models with leader (s). We will conclude, 
in section 14.31 by revisiting the flocking behavior of the C-S model. 

4.1. Flocking of the proposed model. 

Theorem 4.1. Consider the model for self-organized dynamics ()2.2p and assume that its 
influence function (j) satisfies, 

/•oo 

(4.1a) dv{0) < / 02(r)dr. 

Then, its solution, {(xj(t), Vj(t))}j, converges to a flock in the sense that (|1.7p holds. In 
particular, there is unconditional flocking if (j)^ has a diverging tail, 

/oo 
(t)^{r)dr = +00. 

Proof. Since 4>{dx) < (pij ^ 1) the alignment coefficients Ojj in ()3.ip are lower-bounded by 

^ 0(|xj -X^l) (j){dx) 

E.'A(lxfc-x.l) - N ■ 
We now set 9 to be this lower-bound of the ciifs, 

o{t) .= 

so that the global active set at that level, A(6{t)), include all agents. Thus, as noted already 
in p.gp . X{9) = N, and the global version of our main theorem 13.51 vields. 

j/x{t) < dv{t) 

j^dv{t) < -a(l)\dx{t))dv{t). 
The result follows from lemma [3T71 with Tp{r) = 4>'^{r). □ 



Remark 4.2. Theorem \4-l\ tells us that the model h2. 2\) admits an asymptotic flocking 
velocity, Vqo 

limvj(i) = Voo. 

In contrast to the C-S model, however, our model does not seem to posses any invariants 
which will enable to relate Vqo to the initial condition, beyond the fact noted in remark 
that belongs to the convex hull ^2(0). We can therefore talk about the emergence in 
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the new model, in the sense that the asymptotic velocity of its flock, Vqo, is encoded in the 
dynamics of the system and not just as an invariant of its initial configuration. Whether 
Vqo can be computed from the initial configuration remains an open question. 

4.2. Flocking with a leader. In this section, we discuss the dynamical systems with (one 
or more) leaders. 



Definition 4.3. Consider the dynamical system (13. ip . An agent "p" is a leader if there 
exists /3 > 0, independent of N , such that: 

(4.2) aip{t) > /34>{\xp - for every i. 

In other words, an agent "p" is viewed as a leader if its influence on aligning all other 
agents "z", is decreasing with distance, but otherwise, is independent of the number of 
agents, N. We illustrate this definition, see figured! with the following dynamical system: 
a leader "p" moves with a constant velocity and influences the rest of the agents with a 
non- vanishing amplitude < /3 < 1, 

(4.3a) ^ " ^ = " "^vi^j - Vi), 

where 

/3</>(|xp - Xj|), j =p, 

(4.3b) apj = 0, a. 




Figure 4. The agent p (herder) is a leader in the sense of definition (14. 3p . 
He influences every other agents (sheep) more than a certain quantity 
/3(/>(|xp - Xj|). 

We note that there could be one or more leaders. The presence of leader(s) in the 
dynamical system (j3.ip is of course typical to asymmetric systems. We use the approach 
outlined above to prove that the existence of one (or more) leaders, enforces flocking. 

Theorem 4.4. Let {xj(t), Vj(i)} the solution of the dynamical system (|3.ip and assume 
it has one or more leaders in the sense that (j4.2p holds. Then {xj(t), Vj(t)} admits a 
conditional and respectively, unconditional flocking provided (|4.1ap and respectively, ()4.1bp 
hold. 
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Proof. We set 9 = (3(f)[dx{t))- Then the leader "p" belongs to all active sets, Aj(0), and 
in particular, "p" belongs to the global active set A(^). Thus, \{9) > 1. The inequalities 
([A]) yield 

^dxit) < dv{t) 

j^dv{t) < -af<t>\dx{t))dv{t). 

We now apply lemma [22] with il^if) = 4''^{t) to conclude. □ 

Remark 4.5. If the leader p is not influenced by the other agents, then one deduces that the 
asymptotic velocity of the flock Vqo will be the velocity of the leader Vp. But we emphasize 
that in the general case of having more than one leader the asymptotic velocity of the flock 
emerges through the dynamics of (j3.ip . and as with the model (|2.2p . it may not be encoded 
solely in the initial configuration. 

4.3. Flocking of the C-S model revisited. We close this section by showing how the 
flocking behavior of the C-S model (jl.ip can be studied using the framework outlined above. 
By our assumption, the scaling of the influence function i^(-) < 1, we have 



^^</.(|x,-x,|)<l. 



Hence, we can recast the C-S model (jl.lap in the form ()3.2 

1 

N 

- = rv 

dt 



(4.4) ^ = f^^aijivj -Vi 



^(/>(|xi -XjI), jT^i, 



In this case, aij > (j){dx{t)) /N for j ^ i. Moreover, the same lower-bound applies for j = i, 
because of the normalization (f) < 1: 

Therefore, if we now set 9 to be this lower-bound of the a^j's, 

0{t) .= 

then Ap{9{t)), and consequently, A{9), include all agents, X{9) = N, consult ()3.9p . Theorem 
13.51 vields. 

(4.5a) j^dxit) < dv{t) 

(4.5b) j^dv{t) < -a(t>\dx{t))dv{t). 

Now, apply lemma [3^71 with V'(^) = '^^('^) to conclude the following. 
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Corollary 4.6. Consider the C-S model (jl.ip with an influence function, cj), that has a 
non square-integrahle tail, (jl.Sp . Then the C-S solution, {(xj(t), Vj(t))}j, converges, uncon- 
ditionally, to a flock in the sense that ()1.7p holds. In particular, since the total momentum 
is conserved, (jl.3ap . 

v,(t)*^v, v:=l5;v,(0). 

i 

Comparing the quadratic divergence (I4.1bp vs. the sharp condition for C-S flocking, 
(jl.Sp . we observe that the unconditional C-S flocking we derive in this case requires a 
more stringent condition of the influence function. This is due to the fact that the proposed 
approach for analyzing flocking is more versatile, being independent whether the underlying 
model is symmetric or not. 



5. From particle to mesoscopic description 

We would like to study the model ()2.2p when the number of particles N becomes large. 
With this aim, it is more convenient to study the kinetic equation associated with the 
dynamical system (12. 2p . The purpose of the section is precisely to derive formally such 
equation. 

We introduce the so-called empirical distribution [24j of particles /^(t,x, v), 

1 ^ 

(5.1) /'^(t,x,v) := _^<5,^(,)®5^^(,), 

i=l 

where (5x<^^y is the usual Dirac mass on the phase space x W^. Integrating the empirical 
distribution in the velocity variable v gives the density distribution of particles p^(t,x) 
in space, 

1 ^ 

(5.2) p^(t,x) = -J]5.^(,)(x). 

1=1 

Using the distributions and , the particle system (j2.2p reads, 

(5.3a) ^ = v„ 

dvj _ /y w (t>{\y-^i\) (w-Vj) /^(y, w) dydw 
^ ^ ^ ~dF - /y<^(|y-x.|)p^(y)dy • 

Therefore, we can easily check that the empirical distribution satisfies (weakly) the 
Liouville equation, 

(5.4a) + V . Vx/ + Vv • {F[f] f) = 0, 

where the vector field F[f] and the total mass p are given by, 

f (/)(|y — x|) (w — v) f (y, w) dydw r 
(5.4b) F[/](x,v) := a ^^'-^"; ' ^f' , p(y) = / /(y,w)dw. 

Jy'A(|y-x|)^(y)dy 

To study the limit as the number of particles approaches infinity, we first assume that 
the initial condition /q^(x,v) converges to a smooth function /o(x, v) as N ^ +oo. Then 



16 



SEBASTIEN MOTSCH AND EITAN TADMOR 



(5.5) 



it is natural to expect that /^(t,x, v) convergences to the solution /(t,x, v) of the kinetic 
equation, 

' ^tf + ^r■ Vx/ + Vv • {F[f] f) = 0, 
ft=o = /o 

However, the passage from the discrete system (|5.3p to the kinetic formulation ()5.4p is more 
delicate than in the argument for the C-S model [14^ I15j: here, the vector field F[f] may 
not posses enough Lipschitz regularity due to the normalizing factor at the denominator of 
(|5.4bj) . But since this question does not play a central in the scope of this paper, we leave 
the study of existence and uniqueness of solution of the kinetic equation (j5.4p for a future 
work, and we turn our focus to the hydrodynaniic model. 

6. Hydrodynamics of the proposed model and its flocking behavior 



Having the kinetic description associated with the particle dynamics (j2.2p , we can derive 
the macroscopic limit of the dynamics [111 \T0\ I15j . We also extend our method developed 
in section [3] to prove the flocking behavior of the model in the macroscopic case. To this 
end we extend the notion of active sets from the discrete setup the continuum, and the 
corresponding key algebraic lemma 13.11 for skew-symmetric integral operators. 



6.1. Macroscopic system. To derive the macroscopic model of the particle system ()2.2p . 
we just integrate the kinetic equation ()5.4ap in the phase space. With this aim, we first 
define the macroscopic velocity u and the pressure term P, 

p(t,x)u(t,x) = / v/(t,x,v)(iv , P(i,x)= / (v - u) ® (v - u)/(t,x,v)(iv, 

J V J v 

where p is the spatial density defined previously (I5.4bp . Then integrating the kinetic equa- 
tion (I5.4ap against the first moments (1, v) yields the system (see also |15j). 

(6.1a) Otp + Vx • (pu) = 

(6.1b) 94(/9u) + Vx- (pu0u + P) = S(u), 

where the source term S(u) is given by, (recall the notation of ()1.9bp . (w) = cf)* (wp)), 



'(|y-x|)p(x)p(y)(u(y) -u(x)) / / \( \ 

(6.1c) S(u)(x) = a ^ = ap(x) ( ^^^ff^ - u(x) 



0(|y-x|)/?(y)dy 



(l)(x) 



y 



The system (j6.ip is not closed since the equation for pu ()6.1bp does depend on the third 
moment of / which is encoded in the pressure term P. In order to close the system, we 
neglect the pressure, setting P = (in other words, we assume a monophase distribution, 
/(t,x, v) = p{t,x) 5u(t,x)(v)). Under this assumption, (|6.ip is reduced to the closed system 

(USD, 



(6.2a) dtp + Vx • (pu) = 

(6.2b) dtipu) + Vx • (pu (g)u) = S(u). 

We want to study the flocking behavior of general systems of the form (consult figure [5]), 

(6.3a) dtp + Vx • (pu) = 

(6.3b) 9tu-F(u- Vx)u = a(u-u). 
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The expression on the right reflects the tendency of agents with velocity u to relax to the 
local average velocity, u(x), dictated by the influence function a(x, y), 

(6.3c) u(x) = / a(x,y)p(y)u(y)dy, / a(x, y)p(y) dy = 1. 

J y J y 



The class of equations (|6.3p includes, in particular, the hydrodynamic description of our 
self-organized dynamics model, (j6.2p . with 

(6.4) a(.,y)- ^^^^'^^^ 



I«ly-x|Wy)<iy' 



Supp(p) 



Figure 5. The quantity a(x, y) ()6.4p is the relative influence of the particles 
in y on the particles in x. 

We begin with the definition of a flock in the macroscopic case. 

Definition 6.1. Let p{t,x.) > and u(t,x)) be the density and velocity vector field which 
solve (j6.3p . Let Supp{p) denotes the non-vacuum states, 

Supp{p) := {x e M*^ 1/ p{x) ^ 0}, 

and consider the diameters, dx{t) and dy[t), of p and, respectively, m, 

(6.5a) dx{t) := sup{|x-y|, x, y e ^^^^(^(t))}, 

(6.5b) dvit) := sup{|u(t, x) - u(t, y)| , :x.,y e Supp{p{t))}. 

The solution {p, u) converges to a flock if its diameters satisfy, 

(6.6) sup dx{t) < +00 and lim dv{t) = 0. 

Clearly, in order to have a flock, the initial density, po, needs to be compactly supported. 
Furthermore, we also impose that the initial velocity, uq, has a compact support, assuming: 

(6.7) dx{0) < +00 and dv{0) < +oo. 

In the following, we assume there exists a smooth solution (p, u) of the system (j6.3p . 

Hypothesis: Consider the system ()6.3p subject to compactly supported initial data, (po, uq), 
(j6.7p . We assume that it admits a unique smooth solution {p{t),u{t)) for all t>0. 
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6.2. Active sets at the macroscopic scale. To prove that the solution (p, u) converges 
to a flock, we need to show that the convex huh in velocity space, 

Q{t) := Conv{u(t,x) j x € Supp{p{t, •))}, 

shrinks to a single point, as its diameter, dv(t), converges to zero. To this end, we employ 
the notion of active sets which is extended to the present context of macroscopic framework. 
We begin by revisiting our deflnition of active set using the influence function a(x, y) (j6.4p . 

Definition 6.2. Fix 9 > 0. For every x in the support of p, we define the active set, Ax(^), 
as 

(6.8) Ax(0) = {y € Suppip) I a(x, y) > 8}. 

The global active set A{6) is the intersection of all the active set Ay^{9): 

(6.9) A{e) = Pi Ax(6') = {y € Supp{p) \ a(x,y) > 9 for all x in Supp{p)] . 

As before, we let Xi{9) denote the density of agents in the corresponding active set; thus 

(6.10) Ax(^):=/ p{y)dy, Xi9) = [ p{y)dy. 

JAx(e) JA(6») 

We would like to extend the key lemma 13.11 from the discrete case of agents to the 
macroscopic case of the continuum. This is formulated in terms of the maximal action of 
integral operators which involve antisymmetric kernels, A;(x, y). 

Lemma 6.3. Let p G L^(M.'^) be a positive function and let k be a bounded antisymmetric 
kernel, |/c(x,y)| < M and /c(x, y) = — A;(y,x). Fix two positive functions, u and w in 
with a total mass U and W, 

(6.11) ^= / ^(x)p(x)dx and W= u;(x)p(x)dx. 

Then, for every positive number 9, we have: 



.12) 



< MUW {I - X^{9)9^). 



k{x, y)u{x)w{y) p(x)p(y) dxdy 
Here, X{9) is the density of active agents at level 9 for u and iv, 

X{9) = I /3(x)dx, Au,w{d) ■■= {x € Supp{p) I n(x) > 9U and w{x) > 9W}. 

Proof. To simplify, we denote 5 := ^ A;(x, y) ti(x)u;(y) p(x)p(y) dxdy. The anti-symmetry 
of k enables us to rewrite, 

\ I H^^y) u{^)w{y) - u{y)w{x) p(x)p(y) dxdy. 

The bound on k and the identity \a — b\ = a + b — 2 min(a, b) yields, 

|5| < / M u(x)w(y) + ii(y)'u;(x) - 2 min ('u(x)i(;(y), ii(y)'u;(x)) /3(x)/)(y) dxdy. 

Using the notations, we obtain, 

\S\<MUW — M / min ('u(x)i(;(y), u(y)t(;(x))p(x)p(y) dxdy. 



A NEW MODEL FOR FLOCKING 19 

We now restrict the domain of integration on the right hand side to (x, y) € x Ag, where 
the lower-bounds of u and w yield, 

|5| < MUW -M I eueWp{x)p{y) dxdy = MUW - MO'^UWX'^ie), 

and (161^ follows. □ 

6.3. Decay of the diameters. The diameters {dx,dv) also satisfy the same inequality 
at the macroscopic level. We only need to adapt the proof using the characteristics of the 
system (|6.3p . 



Proposition 6.4. Let (p, u) the solution of the dynamical system ()6.3p . Fix an arbitrary 6 
and let X(6) be the density of agents on the corresponding global active set A{9) associated 
with this system, ()6.10p . Then, the diameters dx{t) and dv{t) in (|6.5p satisfy, 

(6.13a) ^dxit) < dvit) 

dt 

(6.13b) - -a>^'^i^)^^it)dvit)- 

Proof. We fix our attention on two characteristics X{t) = u{t,X) and Y{t) = u{t,Y), 
subject to initial conditions, X{0) = x and 1^(0) = y for two points (x, y) in the support 
of p(0). Their relative distance satisfy: 

^\Y -Xl"^ = 2(Y-X , u(Y)-u(X)) < 2dxdv. 
dt 

Since this inequality is true for every characteristics, (|6.13ap follows. 

We turn to study the relative distance in velocity phase space: using (|6.3cp we find, 

I |u(y) - u(x)|2 = 2a (u(y) - uix) , u(y) - u(x) - (u(y) - u(x))), 

and hence, 

(6.14) ^ |u(y) -u(X)|2 < 2ady(|u(y) - u(X)| - \u(Y) - u{X)\ 

Using the fact that a{X, •) has a unit p-mass, the difference of averages u{Y) — \i{X) can 
be expressed as 

u{Y)-vl{X) = I a{Y,w)p{w)u{w)dw - I a{X,z)p{z)u{z)dz 

a(y, w)/9(w)u(w) a[X, z)p{z) dwdz 
a{X, z)p(z)u(z) a(Y, w)p{w) dwdz 

a(y, w)a(X, z) [u(v^) — u(z)] p{w)p{z) dwdz. 

We now appeal to the maximal action lemma, 16.31 with anti-symmetric kernel, A;(w, z) = 
u(w) — u(z), and the positive functions n(w) = a(y, w) and w{z) = a{X,z): since 

|u(y)-u(z)| <dy , U= [ a{Y,y)p{y)dy = 1 and = / a(X, y)p(y) dy = 1, 
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we obtain, 

|u(y)-u(x)| <dv{i-\\e)e^), 

Inserted into (j6.14p . we end up with 

^ |u(y) - u(X)|2 < 2adv[dv{l - \^{d)e^) - |u(y) - n{X)\ 

Finally, since the support of p is compact, we can take the two characteristics Y{t) and 
X{t) such that at time t we have: dv{t) = \u{Y) — u{X)\, and the last inequality yields 
(l6libll . 

d_ 
It 



■dl{t) <2adv{t)(dv{t){l- X'^{9)e^) - dv{t)^ = -2a\^{d)e^dl{t). 

□ 



6.4. Flocking in the hydrodynamic limit. Since the diameters dx and dy satisfy the 
same system of inequalities at the macroscopic level, ()6.13p . as in the particle level, (|3.10p . 
we immediately deduce that theorem 14.11 is still valid for the macroscopic system (j6.3p . 



Theorem 6.5. Let (/>, u) the solution of the system ()6.3p . If the influence kernel, (j), decays 
sufficiently slow, 



oo 

2/ 



(6.15) / (t)^{r) dr = +00, 

Jo 

then (p, u) converges to a flock in the sense of definition \6.1{ 

Proof. For every x and y in the support of p{t, •), we have, 

/ N 0(|y-x|) (l){dx) 4>{dx) _ f u \^ - f f \A 

Jy0(|y -x|)p(y)dy Jy 0(O)p(y) dy p Jy Jy 

Thus, if we take 9{t) = (p{dx{t))/'p, every point y in the support of p{t, •) belongs to the 
global active set A{6). Therefore, for this choice of 6, we have. 



\(9) = [ p{x) dx= [ p(x) dx = p. 

JA(6») Jsupp(p) 



We deduce that, 



^dx{t) < dv{t) 
dt 

j^dv{t) < -a(l)\dxit))dv{t). 



To conclude, we apply lemma \37l\ with ip{r) = (p'^{r). □ 

7. Conclusion 

There is a large number of models for self-organized dynamics [H EJ HJ O [71 [I3l [I2l [171 [El 
[20l [22l [211 [26l [25l [27] . In this paper we studied a general class of models for self-organized 
dynamics which take the form ()3.ip . 



dXj dVj ■r-^ 

^ = Vi, ^ = « (vj - Vi), > 0, 2^ ay = 1. 
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We focused our attention on the popular Cucker-Smale model, [SI [9]. Its dynamics is gov- 
erned by symmetric interactions, aij = (pij/N, involving a decreasing influence function 
(pij := (/'(Ixj — Xj\). Here we introduced an improved model where the interactions be- 
tween agents is governed by the relative distances, aij = (t>ij / Ylk^iki which are no longer 
symmetric. To study the flocking behavior of such asymmetric dynamics, we based our 
analysis on the amount of influence that agents exert on each other. Using the so-called 
active sets, we were able to find explicit criteria for the unconditional emergence of a flock. 
In particular, we derived a sufficient condition for flocking of our proposed model: flocking 
occurs independent of the initial configuration, when the interaction function (j) decays suf- 
ficiently slow so that its tail is not square integrable, (jl.Sp . Similar results holds for models 
with one or more leaders. This is only slightly more restrictive than the characterization 
of unconditional flocking in the symmetric case, which requires a non-integrable tail of 
(113]). _ 

In either case, these requirements exclude compactly supported i;^>'s: unconditional flock- 
ing is still restricted by the requirement that each agent is influenced by everyone else. A 
more realistic requirement is to assume that (j) is rapidly decaying or that the influence 
function is cut-off at a finite distance. Here, there are two possible scenarios: (i) condi- 
tional flocking, namely, flocking occurs if dy(0) and dx{^) are not too large relative to the 
rapid decay of (/>^, dy(0) < f^^^^ (j)'^{r)dr; (ii) a remaining main challenge is to analyze 
the emergence of flocking in the general case of compactly supported interaction function 
4>. Clearly, this will have to take into account the connectivity of the underlying graph Q, 
()3.8|) . We expect that the notion of active sets will be particularly relevant in this context 
of compactly supported (/>'s. The main difficulty is counting the number of "connected" 
agents in the corresponding active sets. As a prototypical example for the difficulties which 
arise with both — asymmetric models and compactly supported interactions, we consider 
self-organized dynamics which involves vision, where each agent has a cone of vision, 

(7.16) '^=Vi, ^ = K{iJi,Xj - Xj) 4>ij{vj - Vj). 

i=i 

Here, K{u!i,Xj — Xj) determines whether the agent "i", heading in direction Wj := Vj/|vj|, 
"sees" the agent "j": 



1, ifu;,- ,"'^' ^\ >l>-l, 



0, otherwise. 



with 7 being the radius of the cone of vision (see figure [6|) . The (pijS determine the pairwise 
alignment within the cone of vision, and can be modeled either after C-S (jl.ip . 



~ -Xil), Ni := #{j I K{u}i,Xj -Xj) = 1}, 



1 

or after our proposed model for alignment, (j2.2p 

_ (/"(IXj-Xil) 
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In either case, the resulting model (|7.16p reads, 

= a ^ aij{-Vj - Vj), aij = K{uji,^j - Xj) (t)ij, 
i=i 

and it lacks symmetry, aij 7^ aji. The loss of symmetry in this example reflects possible 
configurations in which agent "i" "sees" agent "j" but not the other way around. This 
example demonstrates a main difficulty in the flocking analysis of local influence functions, 
namely, counting the number of active agents aij > 6 inside the cone of vision. We leave 
the flocking analysis of this example to a future work. 




Figure 6. Adding a cone of vision in the C-S model ()7.16p breaks down the 
symmetry of the interaction. Here, the agent "i" does not "see" the agent 
"j" whereas the agent "j" sees the agent "i" . 
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